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Abstract 

Helmholtz’ decomposition theorem for vector fields is presented usually with too strong restric¬ 
tions on the fields. Based on the work of Blumenthal of 1905 it is shown that the decomposition 
of vector fields is not only possible for asymptotically weakly decreasing vector fields, but even 
for vector fields, which asymptotically increase sublinearly. Use is made of a regularization of the 
Green’s function and the mathematics of the proof is formulated as simply as possible. We also 
show a few examples for the decomposition of vector fields including the electric dipole radiation. 
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I. INTRODUCTION 


According to the Helmholtz’ theorem one can decompose a given vector held v{x) into a 
sum of two vector helds vi{x) and Vt{x) where vi is irrotational (curl-free) and h) solenoidal 
(divergence-free), if the vector held fulhlls certain conditions on continuity and asymptotic 
decrease (r —)■ cxo). Here x is the position vector in three-dimensional space and r = |a;| its 
absolute value. 

The two parts of the vector held can be expressed as gradient of a scalar potential 
and curl of a vector potential, respectively. Concerning the validity, the uniqueness of the 
decomposition and the existence of the respective potentials one hnds diherent conditions. 

The fundamental theorem for vector helds is historically based on Helmholtz’ work on 
vorticesi*^ and therefore also known as Helmholtz’ decomposition theorem. For hydrody¬ 
namics this theorem is of particular relevance, since the huid helds of the decomposition 
have the physical properties of freedom of vorticity and incompressibility, which for each 
held makes the analysis simpler. Especially for the visualization of vector helds the decom¬ 
position theorem is of importance^. 

Foppl^ introduced the decomposition theorem into electrodynamics. He assumed a hnite 
extension of the sources and vortices and therefore assumed a behavior for the corresponding 
vector held of the form |F| ~ 1/r^ for |x| = r —)■ oo. However, his proof allows less restrictive 
conditions, namely an asymptotic decay of the held only somewhat stronger than 1/r. The 
decomposition theorem can be found in one of these formulations in most textbooks or 
lecture notes on electrodynamics. 

Already in 1905 Otto Blumenthal^ proved, that any vector held, that goes to zero asymp¬ 
totically can be decomposed in a curl-free and a divergence-free part (weak version). His 
formulation reads as follows^: 

” Let F be a vector, which is in addition to arbitrary many derivatives everywhere hnite 
and continuous and vanishes at inhnity with its derivatives; then one can decompose this 
vector always into two vectors, a curl-free F and a divergence-free Vt, such that 


V = vi{x) + vt{x). 


n 


The vectors vi and Vt diverge asymptotically weaker than Inr. 

In addition one has the following proposition for uniqueness: vi and Vt are unique up to an 
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additive constant vector, because of the given properties.” No further specihcation for the 
behavior of the vector held was given. 

This formulation was taken over in its essential statements by Sommerfeld^ in 1944. He 
noted further that the fundamental theorem of vector analysis, as he called it, was already 
proven by Stokes^ in 1849 and in a more complete form by Helmholtz’ paper of 1858. 

The extension to a decay of 1/r and weaker is important for electromagnetic radiation 
but also for a few conhgurations in electro- and magnetostatics. 

Later on it was shown that the conditions of continuity and differentiability can be 
weakened^i^ and that the theorem can be applied to vector helds behaving according to a 
certain power law^^. Based on Blumenthal’s method of regularization of the Green’s function 
Neudert and Wahl^ investigated among other things the asymptotic behavior of a vector 
held V if its sources div v and vortices curl v fulhll some conditions including diherentiability 
and asymptotic decay. 

These developments remained to a large extent unnoticed in the physical literature^ and 
in mathematical physics^. Thus it was necessary to show the validity of the decomposition 
theorem for electromagnetic radiation heldsl^ that decay asymptotically with 1/r . 

First we develop a systematic method, the so called regularization method, which is the 
basis of Blumenthal’s proof, but is not explicated in its generality and its improvement in 
order to be applicable to vector helds, which decay asymptotically with a specihed power 
law. Then we reformulated the decomposition theorem including all potentials for such 
cases. It is shown how the levels of the regularization modihes the validity of the uniqueness 
for the vector helds depending on their behavior at inhnity. The equations necessary to 
construct the decomposed parts are presented. The necessity of a formulation of a proper 
asymptotic condition either for the irrotational or solenoidal part is pointed out. 

In the next section we apply the decomposition theorem to the electromagnetic dipole 
radiation making clear why the traditional theorem is applicable and how the diherent 
quantities in the theorem are related to electrodynamics. 

Then we formulate and proof the extension of the theorem and uniqueness up to a constant 
vector for sublinearly diverging vector helds. A mathematical example without physical 
background is given in an Appendix. Such a case might be of interest for physics, if one has 
sources (circulations) that remain hnite even at inhnity. In the conclusion we summarize 
the results. 
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II. REGULARIZATION METHOD 


The solution 0o(^) of the Poisson equation 

A(j)o{x) = -47rp(f) 


( 1 ) 


with the source density p{x) is found by introducing it’s Green’s function 


Go(x,f') 

(j)o{x) 


1 


\x' — x\ 

J d^x' p{x')Gq{x,x') . 


( 2 ) 

( 3 ) 


If the solution exists in the whole domain of the integral has to be hnite. This is 
guaranteed by a sufficient decay of the integrand, either by a sufficient strong decay of the 
source density and/or by a sufficient decrease of the Green’s function. 

In his work on the Helmholtz’ theorem^ Blumenthal presented a method to make this 
solution hnite (regularizing the solution) by changing the Green’s function of the Poisson 
equation, without changing the Poisson equation (that means without changing the source 
density). Thus one can prove the existence of the potential for cases where the source density 
is less strong decreasing. From this method it becomes clear how a systematic extension of 
the decomposition theorem is possible. 

Introducing an arbitrary point xq (apart from the condition that p{xq) is hnite at this 
point; regularization point) and noting that Go{x,x') = Go{x — xo,x' — xq), we expand Go 
in a power series in T — To 


Go(x,x') = 


1 


+ 


(x-Xo) ■ (x' -Xo) 


+o(-G^ 

V \x' — Xol 


( 4 ) 


\X' — Xo\ \X' — Xol 

A stronger decrease for large \x'\ of the Green’s function is now reached by subtraction of 
the corresponding expansion terms. We get the following set of stronger decreasing Green’s 
functions 


Gi{x — Xo, ^ — Xo) = Go{x, x') — 


\X' — Xq\ 


G 2 {x — To, T' — Xo) = Gi(x — Xo, x' — To) 


(T - To) • (T' - To) 


\x' — ToP 


( 5 ) 

( 6 ) 


The asymptotic decrease of these modihed Green’s functions is as l/r'i+®. For i <2 the 
subtracted terms do not change the source density 


AGi{x — Xo, x' — To) = —4:Tt5{x' — x) for 0 < i < 2 . (7) 
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But they allow to extend the range of the validity for which the existence of the potential 
(and the decomposition) can be proven 


0j(a;) = J d^x’p{x') Gi{x — Xq,x'— Xq) and A0j(a;) = —47rp(x) for i<2. (8) 

The solutions 0i(a;) differ only by a (divergence- and curl-free) solution of the Laplace equa¬ 
tion, i.e. (j)o{x) differs from 0i(x) by a constant value and from (j)2{x) by a linear function, 
both depending on xq. 

Trying to extend the range of validity even further one may subtract the next (third) 
term in the expansion (0]) from G2 and obtains 

1 ^ 2 1 

G3{x - Xo,x' - Xo) = G 2 {x - Xo,x' - Xo) - -Ux - Xo) -V) y—-(9) 

/ \x' — Xq\ 

But now G3 fulhlls the Poisson equation 

AG3{x — Xq, x' — Xq) = —A'k\ 5 {x' — x) — 5 {df — Tq)] (10) 

from which it follows, that G3 leads to a solution of a modihed Poisson equation 

A03(x) = -47r[p(f) - p(fo)] • (11) 

Thus the method described here is not suitable for Green’s functions Gi with i > 2. This 
means (as we will see later) that vector helds which increase linearly or even stronger will 
not be decomposed by the regularization method described here. 

Nevertheless one should remark that one can solve the Poisson equation even with G3 if 
one subtracts the solution for the inhomogeneity p(xo) 

03(f) = j dV p{x')G3{x - fo, x' - Xq) + (12) 

We refer to this solution in section IV A 41 
The relations 

VGi+i(f — fo, f'— fo) = —V'Gi(f — fo, f'— fo) for i<2 (13) 

can be derived from ([5]), (|6]) and (Ej). They are used a few times, mainly to compute the 
vector helds vi and Vt and to establish relations between them. 

In the following we will restrict ourselves to the regularization point fo = 0, because the 
Green’s functions are simpler without loss of generality. In this case the potential is hxed to 
0(f = 0) = 0. We will keep this choice in the remaining part of the paper as far as possible. 


5 






III. THE FUNDAMENTAL THEOREM OF VECTOR ANALYSIS 

As already noticed, the formulation of the fundamental theorem rests in its form today 
on the work of Blumenthal. However there are several reasons not to take the formulations 
of Blumenthal resp. Sommerfeld literally. For instance the uniqueness of the decomposition 
into the helds of the sources and vortices, was only shown up to a constant vector. We 
will formulate the conditions in such a form, that a strict uniqueness of the decomposition 
is given. Furthermore in the proof, which will be given, the potentials by which the de¬ 
composed helds are calculated, are part of the theorem (strong version). It is common in 
electrodynamics to calculate the physical helds via the introduction of potentials. Moreover 
since the proof of Blumenthal is somewhat complex and lengthy it is not found in detail 
in textbooks. Therefore a shorter and more compact proof seems to be useful. Thus we 
formulate the theorem in the following way: 

Let v{x) be an everywhere continuous differentiable vector Geld with the asymptotic behavior 
lim v{r)r^ < oo, where e > 0, then the decomposition 

r—>-oo 

v(x) = vi + Vt = —V(/)(x) -I- V X A(x) (14) 

is unique with 

A{x) = ^ / dV (V' X - A)- (16) 

Att J ^\x'-x\ r'' 


Remarks: 

• Curl- and divergence-free helds Vh can be added to vi if they are subtracted from 
Vt without ahecting the boundary conditions of v. Such harmonic vector helds are 
suppressed if one explicitely demands that h) and/or Vt vanishes asymptotically and 
establish a strict uniqueness of the decomposition. 

• Usually the potentials (f){x) and A{x) are dehned with the Green’s function Go (E])- If 
they are hnite, then there is no need for Gi (Ej). However if the vector held v decays 
asymptotically as 1/r or weaker, one generally has to use the Green’s function Gi as 
shown in flT^ and ffTbjl in order to avoid divergences in the potentials (j){x) and A{x). 
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• As already mentioned in section m the potentials are fixed to the values 0(0) = 0 and 
A(0) = 0 by the choice of the regularization point xq = 0. This choice does not affect 
the vector fields h) and h). 

• The vector potential A by its definition is purely transversal, V • A = 0 (see proof 
below in App. [A|) . 


• We want to stress the point that the decomposition theorem holds for any vector 
field independent of the type of physical equations the vector field might fulfill. On 
the other hand if one thinks of the electric field or the magnetic field as examples of 
the theorem, due to the Maxwell equations these fields are connected although with 
respect to the decomposition theorem they are independent. 


Let us define the source density p{x) and the vortex density j{x) as 


p{x) = 


V • v{x) 
dvr 


J(^) = 


V X v{x) 
Air 


( 17 ) 


then the decomposition of the corresponding vector field in its irrotational (curl-free) and 
solenoidal (divergence-free) parts leads to the result, that 


V • F(F) 

= Anp^x) 

and 

V X F(x) = 0 

(18) 

V X vt{x) 

= 47rJ(f) 

and 

V • Vt{x) = 0. 

(19) 


1. Proof of the fundamental theorem 

First we show the existence of the scalar potential. If the finiteness of the integral flT5l) 
is proven, one gets the field vi by calculating the gradient of 0. For this it is required that 
the integration and differentiation interchange. Then one can show that V x u) = 0 and 

V • u) = V • F. 

Subsequently one proceeds quite similarly for the vortex field by showing the existence 
of flT^ first, then calculating vt and proving its properties V x F* = V x F and V • Ft = 0. 
Finally we check that the sum vi + Vt = F. 

If we show that the integral flT^ exists and is finite, then the longitudinal part F can 
be determined. We note, that the singularities at x and at zero do not lead to a diverging 
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contribution to the integral, because of the antisymmetry of the integrand around the singu¬ 
larity. More important is the asymptotic behavior of the integral for r' —)■ oo. We integrate 
then over the surface of a larger sphere with radius R. Now we have to take into account the 
regularization term ([S]). Since no assumptions have been made on the asymptotic behavior 
of the sources p = V • u/dvr but only on u, we perform a partial integration. This allows us 
to prove the convergence from the behavior of the vector held v alone. Integrating over the 
volume of the sphere leads to 



The radius R can be chosen in such a way, that the held becomes small. Then the surface 
integral vanishes as l/R'^ and it remains to show convergence of the volume integral. 

In order to achieve this we separate the volume of integration into an inner volume of a 
sphere Sr with radius R^ r and the outer domain r' > R 


0(x) [ dV (u(f') ■ V')Gi(f, f') -1- (j)a{x). (21) 

JSr 

The contribution of the outer domain to the potential has been indicated by 4>a{x)- For an 
estimate of this term one can take the the Taylor expansion of Gi, ([5]), and hnds 


l0a(F)| = 

< 


-1 

Att 

Jr 


'r'>R 




( 22 ) 




ArVr 


eR^ 


Thus the contribution of the outer domain to the potential vanishes as 1/eR'^, and the 
existence of 0(x) has been proved. 

It should be proven that the negative gradient of 0 flT^ represents the curl-free part h) 
of V 


Vi(x) = -V0(f) = ^ j (V' • v{x')) VGi(T, x'). (23) 

Since vi is calculated from a potential curlh) is zero. 

Now it should be shown that vi has the same sources as v 

V • h) = ^ j dV (V' • F(f')) AG'i(f, x') = V-v. (24) 
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Here we used the property ([7]) of the Green’s function. Both vector helds have indeed the 
same sources. In fl2^ one can replace VGi{x,x') by —VGq{x,x') and one obtains the 


longitudinal vector held in a manner that is known from the potential-theory in electro- and 
magnetostatic 



(25) 


The proof for the vector potential goes along the same lines and is shifted to appendix A. 

2. Proof of uniqueness 

Now, the existence of the potentials (ITH]) and (ITB]) has been proven. In the last step of 
the proof, the uniqueness of the decomposition has to be shown. 

We have decomposed the vector held n in a source held vi and a vortex held fy, under 
the boundary condition that the total held |n| vanishes going to inhnity. In order to reach 
uniqueness of the decomposition we demand that |h;| and in consequence also \vt\ vanish 
going to inhnity. 

Assume two diherent decompositions of the vector held v = vi + vt = and consider 

the diherences of the source and vortex helds. Then the diherence of the longitudinal vector 
held Vd = vi — ifi is Si divergence- and curl-free held that vanishes at inhnity. It can be 
derived from a scalar potential (fd, which fulhlls the Laplace-equation (harmonic function). 
The only solution for the potential allowed would be a constant (see also the argumentation 
for fH5l) and (061)). Thus the vector held Vd has to be zero and the decomposition is unique. 

IV. APPLICATIONS IN ELECTRODYNAMICS 

A. Static fields 

In electrodynamics the fundamental theorem of vector analysis is used especially (al¬ 
though not always mentioned) in magneto-statics for magnetic fields in matter—. We con¬ 
sider a permanent magnet, where no volumen current is present. There (in the Gaussian 
system) the magnetization AnM corresponds^ to the vector field v in the decomposition 
theorem, the magnetic field —H to the longitudinal (irrotational) part vi and the magnetic 
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induction B to the transversal (solenoidal) part Vf. The source is given by pn = —V • M 
and the circulation by jn = x M. 

These fields are related by the material equation AirM = B — H. This exactly corresponds 
to the decomposition theorem. For such a case the sources and vortices are near the surface 
of the magnetic body, since the magnetization inside is almost constant. In any case the 
sources and vortices are localized to a finite region and in consequence the corresponding 
source and vortex field decay asymptotically at least as 1/r^. The total vector field v of the 
magnetization is zero outside the magnetic body. 

A quite similar situation occurs in electrostatics in a medium with spontaneous polar¬ 
ization, where no free charges are present^. There the vector field v corresponds to the 
polarization dvrP, the source field vi to the electrostatic field —E and the vortex field Vt 
to the dielectric displacement field D. Again the decomposition therem corresponds to the 
material equation 47rP = D — E. 

Even in electro- and magnetostatics configurations with slow decreasing fields exist. The 
electric field of an infinite straight wire, which bears an electric charge, decays as 1/P, 
where p is the distance to the wire. If on the other hand the wire carries a current, then 
the magnetic field decays as ~ 1/p. In both cases a regularization is appropriate to get the 
potentials from finite integrals over the sources. 


B. Time dependent fields: the electric field of an oscillating dipole 

Periodically moved charge densities p(x, t) = p{x) of frequency u emit a radiation 
field of the same frequency. For simplicity we use the complex notation understanding 
the physical quantities (charge density, potential, fields) always as the real parts of the 
corresponding complex quantities. The radiation fields factorize in the same way as the 
sources v{x,t) = v(x)e~^^, where in v{x) the dependence on the frequency uj resp. wave 
number k = u/c has been suppressed. A decomposition of the time independent vector field 
v{x) is possible, since the radiation field, or more precisely its long range part, decays as 1/r 
and thus fulfills clearly the conditions of the decomposition theorem. 

If one starts from the assumption that the asymptotic behavior of the field has to be 
stronger than 1/r, additional considerations are needed in order to proof the decomposition 
of the radiation fieldsl^. Radiation fields, which decay asymptotically as 1/r are rarely con- 
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nected with the decomposition theorem. One reason might be that in most of the textbooks 
on electrodynamics the result of Blumenthal’s proof are not mentioned and one gets the 
impression the decomposition theorem can only applied under additional conditions'^ as 
they are found in radiation helds like e*^'’/r. The peculiarity of these cases is, that one does 
not need the regularization term, although one has a held of 0{l/r). 

Strictly speaking the conditions of the theorem are not fulhlled if the vector held has 
singularities due to point sources. This also holds for the radiation helds considered. However 
the integration over the sources in flT^ and ffT6|) remain hnite. The only consequence, in 
cases where a regularization is necessary, is that the regularization point has to be diherent 
from the singular points due to the source. 

The electric radiation held E{x) = v{x) of an oscillating point dipole p(f) = p reads^^ 


vix) =-< X (p X Cr) H—r(l — ikr) 3(p- eAcr — p \ 

r I L \ ) 


(26) 


Cr = x/r is the unit vector in the direction of x and v{x) is the spatial part of the electric 
held. For fc = ca/c = 0, one obtains of course the static dipole held. Let us hrst calculate 
the source and vortex density 


Ph{S) = -ikr)pp{x)= pp{x) pp{x) =-p ■V5{x) (27) 

jh{x) = ^ -(1 - tkr) (e- X p) = -B{x). (28) 

Pp{x) is the localized charge density of the static dipole^, whereas the vortex density is 
extended in the whole domain decreasing for r —)■ oo as the radiation held with 1/r. It can 
be identihed with the spatial part of the magnetic radiation held^ B apart from a factor, 
as expected from Faraday’s law of induction. Surprisingly the wave number dependence in 
Ph{x), which in jnix) comes from the retardation, drops out. This asymmetry has already 
been discussed by Brill and Goodmai>2S. Hence the scalar potential is given by 

= f dV If j = 0c(^)- (29) 

J \x — x\ 

Multiplying by the factor one obtains the quasi-static (acausal) dipole potential 0c(^, i) 
as it it known using the Coulomb gauged. One obtains the spatial part of the scalar potential 
in Coulomb gauged, which the static (acausal) dipole potential. From that it is clear that 
the longitudinal decomposed vector held vi is the quasi-static electric held of a point dipole 


Vi{x) = -V0/^(f) =[-p + 3(p • e;)e;] ■ 


(30) 
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and does not contribnte to the electromagnetic radiation, which is pnre transversal. The 
decomposition is hnally shown by calcnlating the transversal part Vt = V x Ah according 
to the theorem from of the vector potential 






B{x) + X p. 


(31) 


In electrodynamics one never dehnes a vector potential for the electric held, bnt it is known 
from the Ampere-Maxwell-eqnation that the electric held can be calcnlated via the cnrl of 
B. The longitndinal part of vi is removed by the second term of Ah, a qnasistatic vector 
held. Note that this is not the vector potential Ac known from calcnlating the electric and 
magnetic helds in the Conlomb gange 

(1 - } 

= ^{E{x) + V4>c{S)). (32) 


Akr 


Ac{x) = -< ik Cr X {p X Cr) — [p — 3(er • p)e)-] [—h 


l 

kr'2 


Thns with the transverse held 


Akr 


Vt{x) = - < k^Cr X (p X Cr) + -5-(l — ikr) 3(p • 6^)6^ — p 


hw 


)&r 




= E{x) - Vi{x)- 


(33) 


the cansal character of the total electric radiation held v{x) is restoredS^iS^. 

The same decomposition may be done for the magnetic radiation held, which however is 
trivial since the held is only transversal (see fl28|l h The vector potential fnlhlls V ■ Ah = 0 
and V X Ah = B, the same conditions as for the vector potential Ac in the Colonmb gange. 
We obtain indeed Ah{x) = Ac{x)- The vortex density of the magnetic held is apart from 
a factor given by the same expression as the electric held E{x) of the electric dipole 
radiation 

Jh{x) = X B{x)- (34) 

Thns all the helds, the vector potential Ah{x)j the vortex held i?(x) = V x Ah{x) and the 
vortex density JhA) = V x 5(x)/47r decay asymptotically as 1/r. This is a conseqnence 
of retardation. We also note that the last term in flTHll . the regnlarization term, which 
gnaranties the convergence for a weak decrease of the held as 1/r, is not necessary in this 
case. The integrals converge even withont this terml^. This also applies for other helds like 
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v(x) = p/r. On the contrary, for a vector field like v{x) = erjT the regularization term is 
necessary for reaching convergence, but the regularization point xq has to be different from 
zero. Then we get for the potential 

0(x) = Inro — Inr . (35) 

One may be surprised that all calculations for the radiation held could be performed with 
a regularization on a lower level (Go instead of G\ etc) than expected according to the decay 
of the vector held. One reason lies in the symmetries of the sources and circulations (see 
App. [C]). 

V. DIVERGING VECTOR FIELDS 

A. Supplement to the fundamental theorem of vector analysis 

As already mentioned, the fundamental theorem of vector analysis can be applied to 
asymptotically sublinearly diverging vector helds, if one inserts the faster decaying Green’s 
function G 2 dS]) into flT^ and flT6|) for the computation of 0 and A. 

Let v{x) be an everywhere continuous differentiable vector Geld with the asymptotic behavior 
lim n(r)/r^“^ < 00 , with e > 0, then the decomposition 

r—>-oo 

v(x) = vi + Vt + Vc = —V0(x) + V X A(x) + v(xo) (36) 

with 

0(x) = ^ J d^x' (y' ■ v{x'))G2{x — xo,x'— xq) (37) 

A{x) = ^ [ d^x'y'X v{x'))G 2 {x — Xq^ — Xq) (38) 

dvr J 

is unique apart from a constant vector Geld. 

Remarks: 

• The vector potential A by its dehnition is purely transversal, V • A = 0 (sse proof 
below). 

• The regularization of the Green’s function at a point xq is responsible for the hniteness 
of the integrals flTT)) and fl38|) . The point can be chosen arbitrarily. 
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Curl- and divergence-free (harmonic) fields Vh can be added to vi if they are subtracted 
from Vt without affecting v. 


• Harmonic helds with the exception of constant vector helds Vc can be suppressed if one 
demands that vi asymptotically diverge weaker as linearly. This is shown in section 




• If vi and Vt are calculated with (l37|) and fl38|) . then one obtains the value of the constant 
vector Vc = v{xq) depending on the arbitrary regularization point. 

• If the vector held approaches zero slower than any power law or if it diverges logarith¬ 
mically (as it is the case in Blumenthal’s formulation of the theorem), or if it increases 
sublinearly, then the faster converging Green’s function G 2 has to be applied in 0 and 
A. The price one has to pay for this weaker requirements on the vector held v is the 
loss of the rigorous uniqueness of the decomposition. 

• If one uses the regularized Green’s function G 2 for the case where the vector held v 
decreases stronger, one recovers the unique decomposition of the fundamental theorem 
Cl, since all integrals coming from the regularization terms are hnite. 

• In the special case of the theorem where v approaches zero at inhnity weaker as any 
power of 1/r (the case e = 1), then vi and Vt may diverge logarithmically although the 
sum of the two parts decays to zero^. 

1. Proof of the supplementary theorem 

At hrst one has to show the existence of 0 and A, fl37)l and fl38|) . Goncerning the potentials 
0 and A, their integrand has the same asymptotic decay governed by v{x')G 2 {x') ~ 
as 0 and A in the former proof for the fundamental theorem. ^ 2 ( 0 ;, x') ([6]) has compared to 
Gi{x,x') an additional singular term at Xq = 0. We have to prove that the contribution of 
this singularity to 0 (and A) is hnite. For this purpose we integrate over a small sphere of 
radius r/ —)■ 0 around zero = cos f}') 



-1 

47r 




d^x' (V' • v{x')) = —27rp(0) r 


/ 


dr' 


L 


-1 


de'r = o. 
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Now we can be sure that (p{x) and A{x) exist. Starting from ffT^ . we can reformulate all 
equations up to flA6p by replacing Gi by Gj+i. 


Before the decomposed vector helds are computed, one should compare the scalar potentials 
(fT6|) with (l38|) . Because of the use of G 2 in 0 these both potentials differ in linear function 
in X. This applies even to the difference between A and A and has the consequence that vi 
and Vi are indeterminate by a constant vector. Starting with fl37p we build vi as the negative 
gradient of (j) and check if vi has the same sources a v 

vi{x) = —V0(T) = —^ j d^x' (V' • v{x'))VG2{x — xq, x' — xq) (39) 

V • Vi{x) = —^ j d?x' (V' • v{x'))/S.G2{x — Xq, x' — Xq) = V ■ v{x). 

We hnd that everything holds as expected. Now we rewrite vi by using flT^ and perfoming 
a partial integration 



As can be seen from fITOl) one gets vi{xo) = 0. If one compares vi{x) computed with Gq in 
one sees that the additional term of Gi subtracts a (divergent) constant held from the 
hrst term to hold vi{x) hnite. This divergence- and curl-free vector held does not contribute 
to the source density. 


A partial integration in the vector potential (138P yields to (see flAip i 

d^x'v{3f) X'V'G2 {x — xo,x'— Xo) ■ (41) 

As already mentioned, diherA and A by a vector linearly in x. Even A{x) is purely transver¬ 
sal (use (IT^ and compare the result with (IA2p l 

V ■ A(x) = — f dV (T(f') X V') ■ V'GAx - fo, x' - xq) = 0 (42) 

An J '• 

Now the vortex held is calculated from the vector potential by taking its curl and transform 
VG 2 to -V'Gi 

Vt{x) = V X A(x) = ^ f d?x' (V' X v{x')) x V'Gi(x — xq, x' — xq) . (43) 

4vr J 

Analogous in the case of the irrotational vector vi also fH3ll contains in addition to flA3p a 
constant vector held ehecting that ht(xo) = 0. Now it is shown that the vortices of h) are 
the same as for v. Inserting G 2 into flA4p one obtains that V x h) = V x x. 
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Our interest is directed to the sum vi + Vt since both helds have divergence- and curl-free 
constant vectors. They inhibit that v is the sum of vi + Vt as the following calculation shows. 
If one replace in (IA 6 p G 2 by G 3 and takes Xq 7 ^ 0 one obtains 

Vi{x) + Vt{x) = ^ j d^x' /\{y{^)G‘i{x - Xq,x'- Xo)) = v{x) - v{xo). (44) 

Now we identify the constant vector v(xo) with fy in fl36|) . Thus Vc, vi and h) depend all on 
Xq. 


2. On the uniqueness in the case of an increasing vector field 


Let us start from two different solutions of the decomposition vi^t and v'^ with the same 
sources and vortices respectively. Then the difference Vd = fi — v[ is a irrotational solenoidal 
vector held. 

This vector held can be written as the gradient of a scalar potential (fd that fulhlls 
the Laplace equation Afid = 0 in the whole space. Its most general solution in spherical 
coordinates reads: 

00 I 

fidir, = (45) 

l=Q rn=—l 

where aim and film are coefficients which allow the solution to fulhll the boundary conditions 
and Yim are the spherical harmonics. All film vanish because the zero-point is contained 
within the domain and the solution should be regular. We now calculate the radial harmonic 
hux of the vector held 


Vd ■ Cj. Vdr 


dr 




(46) 


/=1 rn=—l 


When the distance r goes to 00 and one notes that Vd 




1-e 


with e > 0 , then all coefficients 


with / — 1 > 1 — e have to vanish, otherwise this would lead to a stronger divergence of Vdr- 
Thus only the terms with I = 0 und I = 1 remain. Therefore the solution reads: 

ttoo 

m=—l 


(j)d{r, fi, (f) = aooYoo -h ^ aimYimi'd, r 


\/4^ 


w ■ X . 


(47) 


and we obtain 


Vd{x) =-Vfid = w . (48) 

Hence the held Vd = w is unique up to a constant vector. The choice of the regularization 
point inhuences the constant vector of h) only. 
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3. Comment on the application of the supplementary theorem 


We are not aware of a an analytically calculable physical example in this case, however 
for numerical calculations the knowledge of the validity of such a theorem is important. 
There are cases where one does not know always the exact asymptotic behavior of a vector 
held. We have already seen that in the case of electromagnetic radiation because of the 
(symmetry) properties of the held one could decompose the held with the Green function 
Go although generally Gi should be necessary. A few statements on the inhuence of the 
symmetry are made in the appendix O 

Generally we want to point out that the application of the regularization schema is not 
only restricted to the decomposition of a vector held, but may also be important for other 
problems where the solution of a poisson equation (scalar or vectorial) is used. 

If one considers sources that remain hnite at inhnity, the vector held belonging to this 
sources should diverge linearly or stronger. Within our schema one then has tu use the next 
step in the regularization procedure. There the vector held has to be computed with G 2 and 
is determined only up to a vector held linear in x. 

4- Stronger diverging vector fields 

We have already seen that asymptotically strong decaying Green’s functions Gi with i > 3 
cannot be treated in the same manner a those for i <2. We restrict ourselves to i = 3 what 
means that the vector held may increase less than quadratically. In this case the potentials 
are given by (see flT^ i 


03 (^) = ^ 

dPx' p{x')G2,{x — xq,x' — Xq) -|- 

(49) 

M{x) = ^ 

1 (fx j{x)G 3 {x xo,x xo) + 27r ^ r^. 



The last term of 03 and A 3 cancels the contribution to the inhomogeneity caused by G 3 . 
Now, V can again be decomposed in vi and Vt except for a linear vector held that depends 
on the regularization point xq- 
Remarks: 

• Besides our statements to the regularization of weak diverging vector helds high sym¬ 
metric vector helds can be decomposed by the method shown here even if they diverge 
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asymptotic 

region 

Exponent 

Green’s 

function 

vi{r 00 ) 

Unique up to 

V ~ 

0 < e < 1 

G2{x,x') 

Vi ~ 

w 

u ~ 1/r^ 

0 < e < 1 

Gi{x,x') 

Vi{oo) = 0 

w = 0 

V ~ 

0 < e 

Go{x,x') 

vi{oo) = 0 

w = 0 


TABLE I. Different cases of a vector fields v{x), which decay asymptotically to zero or increase 
sublinearily (first and second column), can be decomposed into longitudinal (irrotational, curl- 
free) vi and transversal (solenoidal, divergence-free) parts besides a constant vector field w. In 
order to cover all cases one has to introduce regularized Green’s functions (see (0)) (third column) 
respectively. Also shown is the asymptotic condition on the longitudinal field (fourth column). We 
also indicate the extent of the uniqueness of the decomposition (fifth column). 

stronger than assumed so far. This can be seen from flB4p for i = 0 (and z = 1), where 
the additional terms to Gi in G 2 cancel the contributions to Sq (and Si) for r' > r. 

• Simple examples for this feature are v = ar°‘ and a > — 1 or n = e).r“. 


VI. CONCLUSION 

We have presented a proof of the fundamental theorem of vector analysis (Helmholtz’ 
decomposition theorem) for vector fields decaying weakly and extended to even sublinearly 
diverging vector helds. Contrary to the original proofs we can distinguish between different 
cases. Our results are summarized in Tab. [T1 Note however that not only the decay of the 
vector field is important for introducing a regularization but also its symmetry. This extends 
the presentations of this theorem given usually in textbook on electrodynamics. Especially 
the case of weakly decaying fields has been discussed in the physical literature in the context 
of electromagnetic radiation helds. 

Considering the validity of Helmholtz’ decomposition theorem there is no doubt that the 
theorem can be applied quite generally to electromagnetic helds either static or dynamic. 
This was demonstrated by explicit examples. 


18 












ACKNOWLEDGMENTS 

One of the authors (D. P.) thanks W. Zulehner for helpful discussions. 


Appendix A: Existence of the vector potential and its transverse vector field 


For the vector potential fITBl) we have the problem that we do not to know the asymptotic 
behavior of the vortex density. Therefore we show the existence of the integral ffT6|) in the 
same way as in subsection IIII II and obtain after a partial integration 


A{x) = ^ / d?x' v{x') X V'Gi(T, T'). 


(Al) 


This vector potential A{x) turns out to be a purely tranversal vector potential for which the 
divergence vanishes. In order to show this we use flT^ in flAlj) 

-1 


V ■ A{x) = — / d^x' {v{x') X V') ■ V'Go(T, x') = 0. 


dvr 


(A2) 


The fundamental theorem of vector analysis states that the solenoidal part of v held is 
given by (see flAlD ) 


F^(f) = V X A{x) = ^ / X (V' X F(T'))Gi(f,T'). 


(A3) 


Since h) is calculated from the vector potential divn) is zero. 

Now one has to show that the vortices of Vt are the same as those of the given vector held 
V. Therefore we calculate the curl of flA3p and use the identity (V x V)F= V(V • F) — Av 

V X vt{x) = ^ f dV V X [V X (V' X v{^))Gi{x,x')] 


-1 

dvr 


d^x' AGi{x, S') (y' X v{x')) = V X v{x). (Ad) 

Vt is a pure solenoidal held and its vortices of Vt are the same as those of v. 


In a hnal step it is shown the sum of the irrotational and solenoidal held vi + Vt = v 
equals the given vector held. For this reason we reshape 0(T) (1^ and A{x) (lAlll by 
replacing V'Gi(a;, T') with —VG 2 (T, T') according to flT^ 


ys) 

A{x) 


1 

dvr 

1 

dvr 


d^x' V • v{x') G 2 {x, x') 
d^x'V X v{x') G 2 (T, T'). 


(AS) 
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Now the negative gradient of 0(x) is added to the cnrl of A{x) and the identity (V x V)n = 
V(V • h) — An is used 


vi{x) + Vt{x) = -— / d^x' V(V • v{x')) — V X (V X v{x')) G 2 {x^x') 
dvr ./ L J 


-1 

dvr 


d^x' A(v{x') G2{x, x')) = v(x) 


(A6) 


Appendix B: Example of a diverging vector field 

We want to study the following vector held 

V = a X {cr X d)^ (Bl) 

where a is a constant vector, v diverges as ~ -y/r. It seems to be more convenient to 
determine hrst sources and vortices and then to calculate the helds belonging to these 

P(^) = ^[3a^ + (aJ(^) = • er)e^ x a^. (B2) 

To make the computation of the potentials as simple as possible we use the regularization 
point Xq = 0. Then we get 0 from flTTI) es follows 

If 1 

(j){x) = — / d^x' — i=\3a'^ + {a ■ er'Y]G 2 {x,x') (B3) 

4:71 J 2Vr' 

In the next step we introduce spherical coordinates and we hx the primed coordinate system 
by the unprimed vector x: Cz' = G and perform the integration over the azimuth ip' = 
cos'd'). This leads to a replacement of sin(p' cosip' by zero and cos^ p' = sin^ p' by 1/2 

a ■ Cr' = ttx' \/l — cos (p' + tty/ \/l — sin p' + az'C 
(a • Cr'f = a^^(l - ^'2) + (a • - 1) 

For the calculation of the angular integral of 0 one needs to evaluate this two surface integrals 
{dQ' = di' dp') 

S.{ry) = j-J dn'yG2{x,x') (Bd) 

*S'o(r, r') = (^ - i) 6 '(r - r') (B5) 

S2yr') = ^5o + ^{^d(r-r') + ^d(r'-r)} (B 6 ) 
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Now we get for the scalar potential 


0(f) = J [7a^ - (a • er)^]S'o(r, r') - [a^ - 3(a • er)^]S' 2 (r, r')| 

= ^[-7a'^+ ery]^/r^ (B7) 


The analogous calculation for the vector potential yields 

A{x) = -a X (e7 x d)y/r^. 
9 


(B8) 


In the last step, the calculation of the decomposed vector helds, we get 


vi{x) = —V0(f) = “gl ~ ■ drY]er + 4(a • e7)^e7|\/r 

Vt{x) = V X 74(f) = -|2a X (e), x a) — 4(a • er)a + {a - e).)ej, x {a x Cr 

9 1 


r. 


(B9) 

(BIO) 


Thus we have demonstrated that sublinearly divergent vector helds can be decomposed in its 
irrotational and solenoidal components, both diverging as ~ ^/r. Since f(fo = 0) vanishes, 
it is indeed vi + Vt = v. 


Appendix C: On the influence of symmetry 


In the preceding example we used flB4l) for the computation of 0(f). The ^-functions in 
flB4p cuts all diverging contributions for r' —)■ cxo in flB7p . Similarily a Helmholtz vector po¬ 
tential 74(f) calculated for a circulation with the symmetry j(f) = f{r)p exists for arbitrary 
diverging /(r). 

In the case of the electromagnetic example in Sec. IIVBI the Helmholtz vector potential 
Ah{x) of the electric radiation held (1^ the same effect happens. The required surface 
integral is calculated for a current fl28ll of the form jn = f{r)er x p and reads 


S°Ar,r') = ^ / dn'e 




\x 




x'\ 3 




The integral converges for r' —)■ cxd with 1 /which is stronger as expected from the behavior 
of Go alone, but one gets the same result for the surface integral if one uses Gi instead of 
Gq. This makes clear that the type of regularization necessary depends on the symmetry 
of the vector held v considered and why we did not need Gi in the case of electromagnetic 


21 







radiation. 


^ Helmholtz H 1858 Uber die Integrale der Hydrodynamischen Gleichungen, welche den Wirbel- 
bewegungen Entsprechen J. fiir die reine und angewandte Mathematik vol. 1858, no. 55, pp. 
25-55, and Helmholtz H 1867 On Integrals of the Hydrodynamical Equations, which Express 
Vortex-Motion Philosophical Magazine and J. Science 33 no. 226, pp. 485-512, . 

^ Stokes G. On the dynamical theory of diffraction, Trans. Cambridge Phil. Soc. 9, S. 1, Compl. 
Works vol. II, see p. 10, item 8 

^ Already Stokes^ has 1845 performed such an analysis of the movement of a liquid. This is more 
explictly explained in Lamb’s Hydrodynamics Dover Publ. Reprint 1932 chapterl HI item 30 p. 
31. 

^ See the review by Bhatia H, Norgard G, Pascucci V and Bremer P T 2013 The Helmholtz- 
Hodge Decomposition - A Survey IEEE Transactions on Visualization and Computer Craphics 
19 1386-1404 

® Abraham M 1918 Theorie der Elektrizitdt, vol. 1, 5. Aufl. G.G. Teubner Leipzig, 1. Aufl. 1894 
® Blumenthal O 1905 Uber die Zerlegung unendlicher Vektorfelder Math. Ann. 61 235 
^ See Ref.-^ p. 236 

® Sommerfeld A 1970 Vorlesungen uber Theoretische Physik Bd. H, 6. Aufl. Akademische Verlags- 
gesellschaft Leipzig and 1950 Mechanics of Deformable Bodies: Lectures on Theoretical Physics, 
Vol. 2, Academic Press New York 

® Butzer P and Volkmann L, 2006 Otto Blumenthal (1876 - 1944) in retrospect J. Appr. Theory 
138 1-36 

Tran-Gong Ton 1993 On Helmholtz’s Decomposition Theorem and Poisson’s Equation with an 
Infinite Domain Quarterly of applied mathematics 51 23-35 

Neudert M and Wolf von Wahl 2001 Asymtotic Behavior of the Div-Gurl Problem in Exterior 
Domains Advances in Differential Equations 6 1347-1376 

See for an exception Petrascheck D und Schwabl F 2015 Elektrodynamik Springer Spektrum 
Heidelberg, see p. 221, remark 4 there. 

For an exception see Grofimann S 1981 Mathematischer Einfuhrungskurs fiir Physiker Teubner 
Studienbiicher Physik Stuttgart 


22 



There were several items to clarify for time dependent vector fields, especially the question of 
retardation, its connection to causality and the choice of gauge. 

Miller B P 1984 Interpretations from Helmholtz’ theorem in classical electromagnetism Am. J. 
Phys. 52 (10) 948-950 

Petrascheck D und Schwabl F 2015 Elektrodynamik Springer Spektrum Heidelberg, see p. 221 
Stewart A M 2014 Does the Helmholtz theorem of vector decomposition apply to the wave fields 
of electromagnetic radiation? Phys. Scr. 89 065502 

see^ (8.4.5) and (8.4.6) p. 294. The time dependence is in this case contained in the Fourier 
factor 

The factor e*^’'(l — ikr) is a consequence of the singular structure of pp and can be replaced by 
1 in all calculations. 

Brill O L and Goodman B 1987 Causality in the Coulomb Gauge Am. J. Physics 35 832-837 
The appearance of the quasistatic potentials has led to a discussion on the validity of causality in 
the case of Coulomb gauge. However it was recognized earlier— and confirmed later—, that the 
physical quantities are causal and the decomposition is valid also for time dependent (retarded) 
fields. 

Rohrlich F 2004 The validity of the Helmholtz theorem Am. J. Physics 72 412-413 
Jackson J D 2010 Comment on ’Maxwell equations and redundant gauge degree of freedom’ 
Eur. J. Phys. 31 L79 


23 


